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1.  Introduction  and  Summary. 

Trees  in  a  forest,  stars  in  a  galaxy  and  bacteria  on  a  petri  dish  are  a 
few  examples  wherein  there  is  interest  in  the  spatial  distributions  of  indi¬ 
viduals  over  regions.  In  pursuit  of  better  understanding  of  these  distribu¬ 
tions,  attempts  have  been  made  to  use  probabilistic  and  statistical  tools 
[Eberhardt  (1967),  Pielou  (1969),  Pollard  (1971),  Patil  and  Stiteler  (1974), 

Cox  (1976),  Cox  and  Lewis  (1976)3.  Very  little  effort  has  been  made  to  develop 
coherent  statistical  structures  for  the  locations  of  individuals  in  a  region. 

The  main  purpose  of  this  paper  is  to  propose  such  structures  in  order  to 
avoid  inconsistencies  that  can  be  found  in  the  literature,  to  provide  a  sound 
probabilistic  basis  for  a  variety  of  inference  procedures  that  are  used,  and 
to  permit  the  development  of  better  and  consistent  statistical  methodology  in 
the  area. 

We  show  that,  under  certain  conditions,  the  collection  of  random  variables 
representing  the  numbers  of  individuals  in  subsets  of  a  region  can  be  regarded 
as  a  stochastic  process.  The  models  presented  have  a  stationarity  property.  In 
Section  2  we  present  simple  random  models  which  have,  in  addition  to  stationarity, 
a  property  of  independent  increments .  We  prove  that  simple  models  lead  to  a 
compound  Poisson  process.  In  a  number  of  the  references  cited  above,  Poisson 
and  negative  binomial  distributions  have  been  used  to  describe  the  numbers  of 
individuals  in  subsets  of  a  region.  Our  formulation  shows  that  these  two 
distributions  arise  as  marginal  distributions  in  special  cases  of  simple  random 
models. 

In  Section  3  we  introduce  models  that  are  stationary,  but  do  not 
necessarily  have  the  property  of  independent  increments.  We  show  that  a  negative 

binomial  distribution  proposed  by  Eberhardt  (1967)  to  characterize  the  random 
behavior  of  the  numbers  of  individuals  in  subsets  of  a  region  is  a  marginal 
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distribution  resulting  from  two  different  stochastic  models  that  do  not  have 
independent  increments.  Section  4  is  devoted  to  the  development  of  several 
different  stochastic  models  when  the  region  under  consideration  has  finite 
Lebesgue  measure.  It  is  shown  that  the  binomial  distribution  suggested  by 
Eberhardt  (1967)  can  be  derived  from  all  of  these  models. 

2.  Simple  Random  Models . 

*.  » 

Some  basic  notation  is  needed  to  define  spatial  stochastic  models. 

Throughout  the  paper,  R  will  denote  a  Euclidean  space,  A  will  be  the  Lebesgue 
measure  on  R,  X(A)  will  be  the  number  of  individuals  in  a  subset  Ac  R,  and 
E  will  represent  the  set  of  all  Borel  subsets  of  R  with  finite  Lebesgue 
measures . 


Definition  2.1.  Individuals  are  said  to  be  located  in  R  in  a  simple  random  way  if 
the  following  conditions  hold:: 

(i)  Stationarity:  For  every  set  Ac  I,  the  distribution  of  the  random 
variable  X(A)  depends  only  on  A  (A). 

(ii)  Independent  increments :  If  A^,  . .., ^  (m  >  1)  represent  any  m  disjoint 
sets  in  E,  the  random  variables,  X(A^) ,  ...»  XCA^),  are  independent . 


The  objective  of  this  section  is  to  show  that,  under  conditions  (i)  and  (ii), 
the  collection,  (X(A),  A  e  D,  is  a  stochastic  process  determined  by  a  positive 
number  and  a  sequence  of  numbers  in  [0,  1]  summing  to  unity. 

For  any  two  sets  A^,  in  E  and  nonnegative  integers  ,  k2,  we  have 
clearly 


PCX(A1)  =  1^,  X(A2) 


M  =  l  PCX  (A.  n  A7) 
2  j+t=k.  1  2 

j+**i 


j,  X(Aj  n  £2)  »  r,  X(A1  n  Ay  «  t  J  J 


(2.1) 


Now  let  q  be  an  arbitrary  point  in  R  and  Yq(t),  the  number  of  individuals 
in  a  sphere  centered  at  q  having  Lebesgue  measure  t  a  0,  Yq(0)  =  0.  Then, 
using  (i),  (ii),  we  can  write 


PCXCA^  =  kj,  X(A2)  =  k2]  =  l  PCYqCtp  =  j]P[Yq(t2)  =  r]P[Yq(t3) 

j+l=fe 


iJ,  (2.2) 


where  t^,  t2  and  t^  are  A(A^  n  A2),  \ (A^  n  ff2)  and  X(A^  n  A2)  respectively. 
Without  difficulty,  equations  (2.1)  and  (2.2)  can  be  extended  to  any  finite 
number  of  sets  in  I.  The  general  result  is  stated  in  the  following  theorem. 


Theorem  2.1.  For  individuals  located  in  R  in  a  simple  random  way,  the  joint 
distribution  of  any  finite  collection  of  random  variables  from  (X(A),  A  e  1} 
is  determined  uniquely  by  a  joint  distribution  of  a  finite  collection  of 
random  variables  from  (Yq(t) ,  tiO). 

f 

The  proof  is  omitted.  It  remains  to  show  that,  under  (i)  and  (ii),  the 

collection  (Yq(t),  t  i  0)  is  a  stochastic  process  determined  by  a  positive 

number  and  a  probability  sequence. 

Conditions  (i)  and  (ii)  respectively  insure  that  the  process  Yq(t) 

has  stationary  and  independent  increments.  That  is,  for  real  numbers 

Oit,  st,s  ...it  (ra  a  2) ,  the  random  .variables 
1  z  m 

Yq(ti),  Yq(t2)  -  Yq(tj),  ...,  Yq(tm)  -  ¥q(tm_j)  are  independent  and, moreover, 
Yq(ta)  -  Yq(t4_j)  has  a  distribution  depending  only  on  t^  -  t^,  £*1,  ...,m(t0*0) 
Upon  application  of  a  well  known  characterization  of  such  processes,  it  follows 
[see,  for  example,  Khintchine  (1960),  p.  36]  that  Yq(t)  is  the  compound  Poisson 
process  represented  as 
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where  (N(t),  t  a  0}  is  a  simple  Poisson  process  and  { Z p  i  2  1}  is  a  sequence  of 
i.i.d.  random  variables  on  the  positive  integers,  independent  of  N(t) .  If 
Z^  =  1,  then  Y^(t)  h  N(t),  whereas  if  has  a  logarithmic  distribution,  Y^(t)  has 
a  negative  binomial  distribution,  as  is  well  known.  These  two  marginal  distribu¬ 
tions  of  X(A)  have  been  frequently  used  in  the  literature. 


3.  A  Generalization  of  the  Simple  Random  Models. 

In  this  section  we  present  two  methods  of  developing  a  general  stationary 
stochastic  structure  for  the  random  quantities  {X(A)>.  This  is  done  in  two  steps. 
Firstly,  for  ms  1,  we  present  the  joint  distribution  of  XCAj),  ...,  XCA^ ,  where 
Ap  . ..,  A^  are  any  m  disjoint  sets  of  E.  We  then  show  that  this  set  of  finite 
dimensional  distributions  determines  uniquely  the  stochastic  structure  of  {X(A)>. 

For  this  purpose  we  introduce  some  notation.  Let  {M(t),  t  s  0}  be  a 
stochastic  process  on  the  nonnegative  integers,  stochastically  increasing  in  t, 
that  is,  M(t)  s  M(s)  for  0  s  t  s  s,  and  {Wp  i  i  1},  a  sequence  of  random  variables 
with  positive  integer  values,  independent  of  {M(t)l.  We  assume  that  (M(t),  t  i  0} 
and  {Wp  i  2  1)  depend  on  the  parameters  6^  and  0^  respectively,  ranging  in 
parameter  spaces  P^  and  where  ^  is  a  collection  of  infinite  sequences, 

Pj  and  P2  being  probability  spaces  with  F  and  G  the  respective  probability  measures 


Definition  3.1.  Let  Ap  ...»  A^  be  disjoint  sets  in  E  with  respective  Lebesgue 
measures  tp  ...,  tm,  and  let  kp  ...,  be  nonnegative  integers.  We  define  the 

joint  distribution  of  (X(A^),  i  *  1,  . ..,  m>  as 

m  M(t.) 

PCX(A.)  -  k‘i  r  l,  f  J  n  P[  lx  w.  -  k.]dF(e  )dGCe2). 

11  ’  *1*1  i=i  j-i  y  1  A 


(3.1) 
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Definition  3.2.  With  the  same  notation  as  above,  we  define 


/  /  PC  l  W.  =  k.,  i  =  1,  ....  m]dF(01)dG(92), 

j=r-  ,+1  ^ 

1  2  J  i-l  (3.2) 

i 

where  r.  =  J>  M(t  ),  i  =  1,  ....  m,  r„  =  0. 

1  q=l  q 

The  consistency  of  the  set  of  joint  distributions  of  finite  collections  of 
random  variables  from  {X (A),  A  e  I},  generated  by  disjoint  sets  and  defined  by 
(3.1)  or  (3.2),  follows  from  the  stochastic  structure  imposed  on  (M(t), 
t  >  0,  W.,  i  ^  1).  The  stationarity  condition  (i)  is  clearly  satisfied,  since 
definitions  (3.1)  and  (3.2)  vary  only  with  the  Lebesgue  measures  of  the  re¬ 
spective  sets.  The  generalization  of  equation  (2.1)  to  any  finite  number  of 
sets  in  E,  not  necessarily  disjoint,  provides  a  way  to  extend  definitions  (3.1) 
and  (3.2)  to  any  such  sets.  These  extensions  preserve  the  required  consistency 
aid  stationarity  of  (X(A) ,  A  e  £} .  We  have  proved  the  following: 

Theorem  3.1.  If  the  joint  distributions  of  every  finite  number  of  random 
variables  in  (X(A),  At  I)  are  given  by  the  extensions  of  (3.1)  or  (3.2),  then  J 

j 

{X(A)},  A  €  1}  is  a  stationary  stochastic  process.  j 

It  is  important  to  note  that  (3.1)  and  (3.2)  yield  the  same  marginal  distributions  j 
for  X(A),  A  c  l,  but  they  define  different  processes.  Statistical  methodologies 
based  on  the  two  models  may  be  quite  different  and  this  has  been  ignored  in  the 
literature. 

In  Section  2,  we  proved  that  simple  random  models  depend  on  a  positive  nunfoer, 
say  w,  and  a  probabilistic  sequence,  say  (P^,  i  2  1).  Now  we  show  that  those  model 5 j 
are  particular  cases  of  the  stochastic  processes  presented  by  (3.2).  To  do  so,  let 
{M(t),  t  2  0}  be  a  Poisson  process  with  parameter  y,  and  {W^,  i  2  1}  be  an  i.i.d. 
sequence  of  random  variables  given  by  P(W^  =  i)  =  P^  (i  2  1).  In  addition  we 
take  Pj  and  P^  to  be  sets  containing  only  y  and  the  sequence  (P^>  respectively. 

Now  the  extension  of  (3.2)  reduces  directly  to  a  simple  random  model. 


••rfuw  ixi n  rawujJHD  to  me  . 
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Let  us  assume  that  the  sequence  of  discrete  random  variables  {llL  }  is 
degenerate  at  1,  ^  is  CO,  “),  is  only  the  sequence  (1,  0,  0,  ...)» 
(M(t) }  is  Poisson,  and  F  is  given  by 


F(x)  =  0 


gQ  ^  -By  o-l ,  n 

7TT  /  e  1  y  dy  X  >  0 


a,  B  >  0. 


(3.3) 


x  s  0 


Then  the  marginal  distribution  of  X(A),  A  e  E,  computed  by  (3.1)  or  (3.2),  is 
negative  binomial  with  a  =  a  and  b  =  A  (A).  This  marginal  distribution  for  X(A) 

was  assumed  by  Eberhardt  (1967)  and  by  Patil  and  Stiteler  (1974) . 

4.  Models  for  Finite  Regions  . 

In  the  previous  two  sections,  we  presented  stochastic  models  when  individuals 
were  randomly  located  in  an  infinite  region.  In  this  section,  five  different  models 
are  developed  for  situations  wherein  individuals  are  located  in  a  Borel  set  c_R 
of  finite  Lebesgue  measure. 

Let  Ef  be  the  set  of  all  Borel  sets  contained  in  Rj  and  X^(A),  the  number  of 
individuals  in  set  A,  A  e  E£.  Our  objective  is  to  develop  a  stochastic  structure 
for  the  collection  { X^ (A)  ,  A  e  Ej} .  One  way  to  achieve  this  objective  is  to 
select  a  stochastic  process  developed  in  Section  3  for  the  infinite  region  R,  say 
{X(A),  At  I},  and  condition  it  by  an  event  related  to  the  random  variable  X(R^ . 
Specifically,  let  A^,  ...,  Am  be  sets  in  E£  and  let  k^,  . . . ,  k^  be  nonnegative 
integers.  We  define  the  desired  probability  as 


P[Xf(Ai)  =  kA,  i  -  1,  ...,  ml  =  PCX(A£)  =  k^  i  =  1,  ...,  m[X(Rf)J. 


(4.1) 


The  consistency  and  stationarity  of  the  process  {^(A),  A  e  E^  defined  in  (4.1) 


are  self-evident. 
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For  the  four  remaining  alternative  ways  of  developing  a  stochastic  structure 
for  (Xf(A)},  we  refer  to  the  process  (M(t)},the  sequence  (W^>, the  sets  P^  and  P2, 
and  the  probability  measures  F  and  G, which  were  introduced  in  Section  3.  lve 
define  first  the  joint  distribution  of  XfCAj),  . ..,  xf(Am)>  where  A^,  ...»  Am 
are  disjoint  sets  in  Z £  and  msl,  then  extend  these  definitions  to  any  finite 
collection  of  sets  in  Z £.  Since  the  extension  technique  has  been  used  twice,  it 
is  omitted  from  our  present  discussion.  Now  let  A^,  . ..,  A^  be  disjoint  subsets 
of  Rf,  A(Rf)  =  tQ,  A^)  »  t.,  i  -  1,  ....  m,  and  let  k^,  ...»  km  be  nonnegative 
integers.  We  define  the  desired  probability  in  the  following  four  ways: 

m  M(t.)  M(tn) 

PCX,(A.)  =  k.,  i  =  1,  ....  m3  =  /  /  I  P(  f  W  =  k  |  I  0  W  )dF(e  )dG(02).  (4.2) 
f  1  1  P1xP2  i=l  j=l  3  1  j=l  3  1  L 


PtXfCAp  =  ki,  i  =  1, 


. . . , 


r. 

l 


■  W  Jp(.I 

3=ri-l 


P1xP2 


W. 

>1  3 


M(t0) 

...,  m|  l  W  ) 
j=l  3 

dF(ei)dG(e2)# 

(4.3) 


m  M(t.) 

PCX  (A  )  =  k  ,  i  =  1,  ...,  m3  =  f  /  n  PC  T  W 
*  1  1  P-T2  i=1  3=1  J 


k.|MCto))dF(e1)dG(02). 


(4.4: 


PCXf(Ai)  =  kif  i  -  1, 


. . . , 


r. 

i 


m3  =  / 


P1*P2 


/PC  l 

3=ri-i 


w. 

♦1 3 


where  rQ,  ...»  rm  are  as  in  definition  3.2. 


ki,  i  =  1,  ...»  m|M(t0)) 

dF(91)dG(02). 

(4.5; 


The  stochastic  processes  determined  by  (4. 2) -(4. 5)  have  been  constructed  to  be  both 
consistent  and  stationary.  Definitions  (4.2)  and  (4.3)  or  (4.4)  and  (4.5),  yield 
identical  marginal  distributions  for  Xj(A),  A  e  Ej,  but  define  different  processes. 
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If  we  let  M(t)  be  a  Poisson  process,  {W^},  a  sequence  of  random  variables 
degenerate  at  1,  P ^  and  P ^  singleton  sets,  then  the  marginal  distribution  of 
Xf(A),  A  e  1^,  according  to  each  of  the  definitions  (4.1)  -  (4.5)  is  binomial 
with  X(Rf)  corresponding  to  the  number  of  Bernoulli  trials  and  A  (A)/ A (R-) 
corresponding  to  the  probability  of  success  in  a  single  trial. 

5.  Concluding  Remarks. 

The  selection  of  models  for  stochastic  processes  for  particular  applications 
and  the  development  of  pertinent  statistical  methodologies  have  not  been  addressed 
in  this  paper.  We  have  demonstrated  that  models  proposed  may  be  used  to  yield 
marginal  distributions  assumed  in  the  literature;  they  may  also  be  used  to  avoid 
unwarranted  assunqptions  and  inconsistencies  that  arise.  We  propose  in  subsequent 
work  to  use  the  general  models  of  this  paper  to  devise  improved  statistical 
methodologies  for  problems  involving  the  location  of  individuals  in  a  habitat. 
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